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Abstract 

The hght scalar meson nonet above 1 GeV {i.e. the ag, Kq and /o) are studied within the framework 
of QCD sum rules. In conventional QCD sum rules, the calculated masses of this nonet are degenerate, 
and the mass of Kq is always larger than the oq in contradiction with the observed spectrum. After 
improving the correlation function by including instanton effects, the masses are well separated from 
each other. In particular, our results shows glueball content plays an important role in the underlying 
structure of /o(1500). The decay constants are also discussed. 

PACS numbers: 12.38.Lg, 12.39.Mk, 14.40.Cs, 14.40. Ev 



1 



I. INTRODUCTION 



The SU (3) classification of strongly-interacting particles originally proposed by Gell-Mann 
{l| and Zweig {2^ has been a very successful paradigm in particle physics. It is observed from 
the hadronic spectrum that the results following from this naive quark model (NQM) agree 
better with the heavy meson systems than the light ones. This is understandable because 
the heavy quarks inside a heavy meson are non-relativistic, and hence we can deal with their 
kinematics in the framework of nonrelativistic quantum mechanics as a good approximation. 
However, for light mesons where the light components are relativistic, it is hard to say whether 
the nonrelativistic approximation is applicable. We can see that there is more complexity in 
light mesons than the heavy ones from the observed spectrum. The situation is even worse 
in combination with the proliferation of light scalars and their production in charmless B 
decays. In order to accommodate these light mesons in theories consistent with QCD, models 
beyond the naive quark model have been developed, including glueballs [4], multiquark states 
sl and hybrid states 6]. One hopes that these models can supply some reasonable, or at least 
qualitative, explanation of the observed light mesons. 

The underlying structure of mesons with mass near 1 GeV attract much attention. It has 
widely been suggested that the light scalars below or near IGeV [the isoscalars /o(600), /o(980), 
the the isodoublet Kq{800){ot k) and the isovector ao(980)] form a SU(3) flavor nonet, while 
scalar mesons above IGeV [/n(1370), ao(1450), K*{im) and /o(1500)//o(1710)] form another 
nonet Refs. 0, Q suggest that the light scalar nonet above IGeV can be 

accommodated in the conventional qq model with some gluonic component, while the light 
scalars around IGeV are dominated by qqqq states with some 0"'" qq and glueball states. But 
this interpretation is still far from deciphering the puzzle presented by the light scalars. 

It is obvious that the starting point of all the models mentioned above concentrates on the 
kinematic aspect of the component inside the scalars, i.e., in order to reproduce the spectrum 
in theories consistent with QCD, the complexity of the light scalars is attributed to their 
constituents. Maybe one can refer to this as a kinematics-dependent approach. There is 
another viewpoint we can adopt. We should recognize that in the hadronic region perturbative 



^ Following Rcf. [Sj wc refer to this model as the naive quark model. 
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QCD breaks down, and the nonperturbative aspects of QCD are dominant. It is well-known the 
nonperturbative aspect of QCD is difficult to analyze. The nonzero quark condensate signals 
that the QCD vacuum is nontrivial and has a complex structure. In other word, the dynamics 
in QCD vacuum is very different from the trivial one. It is possible that the complexity 
of the light scalar mesons can be attributed to the enigmatic QCD vacuum, in which the 
particle treated as the excitation of the QCD vacuum from the viewpoint of quantum field 
theory. The nonzero value of QCD vacuum expectation values is one of the main ingredients of 



QCD sum rules [Uj, Q, [3] which deals with the low energy nonperturbative aspects of QCD. 
In conventional QCD rules the physical quantities are expressed by a dominant perturbative 
part and corrections associated with vacuum expectation values of various operators. This 
method works well in many cases, but when we apply this method to the pions, it is difficult to 
obtain reasonable results. This difficulty was solved by introducing instanton contribution into 
the QCD sum rules h^l. Instantons — the nontrivial solution to the Yang-Mills field equation 

n 

[l5l ] — play an important role in solving the puzzle.^. Recent work involving QCD sum rules 



with instanton effects include the electromagnetic pion form factor [17| and glueballs [18|. 
Furthermore, instanton effects within QCD the non-strange sum-rules for scalar currents have 
previously been shown to split the degeneracy between the ao and /o [3].'^ All these works 
pave a new way to resolving the controversy concerning the nature of the light scalars. 

Comparing with the kinematics-dependent approach, we refer to the instanton effects as 
a dynamics-dependent approach, because here one attempts to solve the problem by further 
investigating low energy QCD itself. Keeping these motivations in mind, in this paper we 
investigate the masses of the scalar nonet above IGeV [i.e, /o(1370), ao(1450), -^"0(1430) and 
/o(1500)//o(1710)] from QCD sum rules based on scalar interpolating fields including the cor- 
responding instanton contribution. Because there is no mixing between oq, Kq meson and the 
glueball, these two members are ideally suited to investigate the role of instantons in QCD sum 
rules, and we will analyze them in the naive quark model. The situation is more complicated 
for the /o because it is widely accepted that there is mixing with the isoscalar O"*"^ glueball 
ground state around ISOOMeV 



^ For details on instantons in QCD, see the excellent review by T. Schafer and E. V. Shuryak 16| 
^ We note that these works did not consider the structure of the entire nonet. 
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Because of this mixing, a more consistent analysis should consider the mixing of quark and 
gluonic content in analyzing /o meson. So we will employ a mixed quark-glueball current to 
discuss /o meson if necessary. Specifically, we assign /o(1500) and /o(1710) to be a mixed current 
of quark and gluonic content, while /o(1370) is still assumed to be purely of quark content. As 
will be demonstrated below, the validity of these assignments is upheld by the results of the 
QCD sum-rule analysis. The instanton contributions to the sum-rules are calculated using the 
semi classical approximation with quark zero modes. As a byproduct, the decay constants of 
the states are obtained naturally. 

In Section II we derive the QCD sum rules with scalar interpolating fields in the absence of 
instantons and note the shortcomings associated with the results of this analysis. In Section III 
we present the sum rule including instanton contribution based on scalar current or its mixing 
with gluonic current, and the masses and decay constants of the nonet are calculated. Section 
IV is devoted to our conclusions. 



II. SUM RULES WITHOUT INSTANTONS 

In this section we will discuss QCD sum-rules without instantons and the results following 
it. The starting point is the following correlator defined in terms of scalar interpolating current: 

Ii{q^)=i j rf^a;e^«-^(0|j(x)/(0)|0), (1) 

where j{x) is a scalar composite operator defined as: 

j[x) = qi{x)q2{x). (2) 

compared with the definition in jl^, we have suppressed the renormahzation invariant factor 
(ln(/i/A))~^/'', with /i is the normalization pint and b = (llA'^c — 2?^/)/3. The correlator can be 
expressed in terms of operator product expansion, up to order-a^ perturbative correction and 



Lattice gauge calculations predict a glucball mass of 1400 to 1800 MeV 
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dimension-six, the operator product expansion we get is 

-,2 
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This is the theoretical side of the QCD sum rule from the quark-gluon dynamics point of view. 
On the other hand, the correlator can also be derived phenomenologically: 



1 , imnpi^fs) 

as h subtraction constants . 







s — q^ 



(4) 



The quantity ImnP^(s) obtained by inserting a complete set of quantum states S|?T,)(n| into 
Eq. (1), which reads: 
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where Sq represents the onset of the QCD continuum. The decay constant in Eq. (5) is defined 
as: 

{S\q2qi\Q) = msfs- 

By equating both the theoretical and phenomenological sides, we obtain the total dispersion 
integral: 



1 
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as h substraction constants, 
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(6) 



After Borel transform and subtracting the perturbative continuum contributions, we obtain the 
following sum rule: 
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where scale dependence of decay constant fs is: 



The parameters in Eq. (7) are as follows 24. l25l|: 



a, = 0.517, = 0.012 ±0.006GeV^ 

(uu) = (dd) = -(0.24 ± 0.1)^GeV^ {ss) = (0.8 ± 0.2){uu), 
!Ihi±I!^ = 5MeV, m, = 120MeV, 

{ggUaGu) = {ggdaCd) = O.SGeV'^ {uu) , {gsSaGs) = 0.8{gsUaGu). 

All the values adopted here are given at the scale fi = IGeV. By taking the logarithm of both 
sides of Eq. (7) and applying the differential operator M^d/dM'^ to them, we derive the desired 
mass formula which is free of the decay constant. 

The task now is to find ranges of parameters and the continuum threshold sq such that 
the resulting mass does not depend too much on the value of these parameters. In addition, 
the continuum contribution that is the part of dispersive integral from Sq to oo subtracted from 
both sides of Eq. (7) should not be too large (less than 30% of the total dispersive integral), 
and the contribution of the dimension-six operators is less than 10%. One more requirement is 
the value of the continuum threshold Sq should not stray too much away from the next known 
resonance in that channel 231 ]. 



Before proceeding with our analysis we note that experimentally, except for the /o(1710), 
there is a small mass difference between other members, and the mass difference between qq 
and Kq is even smaller. Thus it is reasonable to deal with them using the same threshold 
and Borel window; we think this criterion also holds true for other multiplets with small mass 
differences considered by other QCD practitioners. Of course one can analyze each member of 
a multiplet with a separate threshold and Borel window, but it is too artificial to be adopted 
because the sum rules are sensitive to the threshold. In considering this, we will analyze this 
nonet within same threshold and Borel window below. 

If one uses the same threshold and Borel window, it is easy to see the mass spectrum of 
the nonet following from Eq. (7) will be similar to the naive quark model; there is a mass 
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FIG. 1: Mass of i^Tg (solid line) and cq meson (dashed line) from sum rule Eq. (7) based on naive quark 
model as function of Borel parameter without instanton. 

degeneracy broken by a tiny difference resulting from the SU(3) flavor symmetry breaking.^ 
Even worse, the mass of Kq with underlying structure sd is always larger than the Oq with 
underlying structure -^{uu — dd). For definiteness, when we select sq = 4.1GeV^ and 
within the range [1.3, 1.6]GeV^,the calculated mass of Kq and Oq is shown in Figure 1. One can 
see the mass of the Kq is above the Oq, which is inverted compared to the experimental results. 
All the results following from Eq. (7) are unsatisfactory. 

To summarize Section II, we conclude that in the conventional QCD sum rule analysis based 
on the naive quark model, one can not separate the this nonet with the same threshold and 
Borel window, the masses following from Eq. (7) are degenerate, and the results for the Kq 
and oq, are in contradiction with experiment. This suggests that important effects have been 
neglected in Eq. (7). 

^ In addition, the degeneracy could also be broken when mass corrections proportional to cis'Ti^j.^-j are taken 
into account. These corrections are negligible compared with instanton elFects which we will consider in the 
next section. 
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III. SUM RULE WITH INCLUSION OF INSTANTON CONTRIBUTIONS 



Basic formula 



It has been known for a long time that the instanton plays an important role in nonper- 
turbative QCD. The starting point on this subject is the solution of classical field equations 
in four dimension Euclidean gauge- field theories given by Belavin et al. [l^. Subseque ntly 
t'Hooft derived the instanton with topological quantum number ri = 1 in Euclidean space 261]: 



^ g [x — xoj + p 
4 p2 



(8) 



where p is instanton size, rjafj.u is the t'Hooft t] symbol, xq is an any point in Euclidean space. 



The density n{p) of instanton with size p in the vacuum can be parameterized as Q, I27I : 

n{p) = ncS{p - pc), 



(9) 



with two parameters and pc, called the average instanton density and size: 

1 



X 10-*GeV*, Pc 



0.6 



(10) 



When we include the instanton contribution in the correlator (1), there is a new term 

2 

d^x e''^-''qio{x)q2oix) 



(11) 



where gio and g2o is the t'Hooft quark zero mode, respectively, m\ and is the effective mass 
correspondingly. Similarly applying dispersive relation to Eq. (11) we can rewrite Eq. (11) as: 



1 , ImH''«'"^^*(s) 
- / ds 



(12) 



/o 5 - r 

After the Borel transformation we get the desired form of the instanton contributions of the 
current with isospin /: 



H 



qq, inst 



.1 n^P 



exp[ 



(13) 
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FIG. 2: Mass of i^Q (solid line) and ao (dashed line) from sum rule Eq. (15) as function of Borel 
parameter include instanton 



and the instanton continuum contribution is: 

where Ko, Ko are the McDonald functions, and Ji , Yi are the Bessel functions. 
When the smoke clears, we get the final result: 

m|/|exp[-m|/M2] = n°P^(M2) - n«™'^°°*(so, M^) 

This is the sum rule we obtained including instanton effects in the correlation function. Similarly 
we can obtain the mass from Eq. (15) with the same manipulation as the previous section. 
As an important parameter in sum rule Eq. (15) it is necessary to discuss the value of the 



effective mass m*'s. In the mean- field approximation 16[ |: 



1 

2\ 2 



m: = m^ = 7rp(-) (iV/\/) ^ = 170MeV , (16) 



9 



The value of N/V is N/V = Ifm ^ phenomenologically. And there is a relation between (uu) 
and {ss): 

{ss) = (0.8 ±0.2) {uu) , (17) 

together with 

m = , (18) 

m* 

we obtain the effective mass of the strange: 

ml = 215t^^MeV. (19) 
In following we take the central value of m*, i.e, m* = 220MeV. 

B. Mass, decay constant of Kq and oq meson 

All the parameters needed in numerical calculation have now been fixed. Firstly we take 
j = ds which corresponds to the Kq meson as our "sample", since in this case there is no 
mixing with glueball, allowing us to examine the effect of instants without the complications 
presented by mixing. . 

Using standard QCD sum-rule methodologies, we obtain the threshold Sq = 3.5GeV^, which 
is just below the next excited state ii'*(1950), and the Borel window is within the range [2.0, 
2.3]GeV^. The calculated mass of is mx* = 1436 ~ 1462MeV. For the oq, assigning 
j = -^{uu — dd) with the same threshold and Borel window, the mass of the Oq is m^g = 
1432 ~ 1472MeV. The results for these two states are shown together in Figure 2. 

We can see from Figure 2 that there is a crossover at the point = 2.16GeV^ corresponding 
to the mass m = 1442MeV, which is very close to the experimental value of rriK* and ao- We 
refer to this value of as the "key point" representing the Borel scale where the mass hierarchy 
between the ao and Kq reverses. The more important aspect observed from figure 2 is that the 
calculated mass for Kq and ao present the right picture : there is a range where the mass of the 
ao is larger than the Kq, a result which cannot be obtained in the sum rule without instanton 
contributions shown in figure 1. In other words, the sum rule including instanton contributions 
can reproduce realistic results which in agreement with the light meson spectrum. 

Having determined the mass, it is straightforward to obtain the decay constant from Eq. (15). 
The decay constants of the Kq and ao are shown in figure 3 and figure 4 respectively. It is obvious 

10 



0.60 



0.55 
0.50 
0.45 

2.00 2.05 2.10 2.15 2.20 2.25 2.30 

M'(GeV') 

FIG. 3: Decay constants of Kq with quark structure sd as function of Borel parameter includes 
instanton. 

that the results are very stable within the Borel window, and from the Figures we find the decay 
constants of Kq and ao: 

Ik (5(1430) 

(IGeV) = SlOMeV, /ao(i450)(lGeV) = 514MeV 

C. Mass, decay constant of /o meson with underlying structure -^{uu + dd) 
Similarly if we set 

j = —^{uu + dd), 
v2 

with isospin / = 0, we immediately get the mass of /o : 

m/o = 1314 ~ 1391MeV, 

At the "key point" i.e, = 2.16GeV^, the mass is m/,j = 1376MeV which is close to the 
experimental value of /o(1370), while the more important result is that the mass of this state 
is no longer degenerate with other states. The decay constant following this mass is: 

//„(i370)(lGeV)= 520MeV. 
11 
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FIG. 4: Decay constants of with quark structure {uu — dd)/V2 as function of Borel parameter 
includes instanton. 

The mass and decay constant of /o are shown in figure 5 and figure 6, respectively. 

Because there are more controversies for the /o meson than the two members discussed above, 
it is useful to mention the model beyond pure quark viewpoint. We also notice the glueball can 
mix with scalar mesons nearby, so it is possible that there is mixing between /o(1370) glueball, 
or in other words, /o(1370) is not in a pure quark state, and could have some glue content; 
an idea that was first introduced in Ref. 



and then generalised in Ref. (29||. The work in 
Ref. [29|] suggested the /o(1710) is dominated by ss content, while /o(1500) and /o(1370) share 
roughly equal amounts glueball(~ 40%) (/o(1500) and /o(1710) will be detailed in the coming 
subsection). Here we see when including the instanton effects, /o(1370) can be accommodated 
naturally in pure quark model. So we conclude that /o(1370) may be a pure quark state based 
when instanton effects are considered. 
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FIG. 5: Mass of /q with quark structure {uu + dd)/^/2 as function of Borel parameter includes 
instanton. 
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FIG. 6: Decay constants of /o with quark structure {uu + dd) / \f2 as function of Borel parameter 
includes instanton. 
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FIG. 7: Mass of /o with quark structure ss as function of Borel parameter without instanton. 

D. Mass, decay constant of /o meson with underlying structure ss 

The /o(1500) may be the most controversial object in this nonet. As we have seen the im- 
portant role of the instanton in giving the mass of the three members in previous subsections. 
Firstly we do some tentative calculations based on pure quark model in hope that these calcu- 
lations shed some light on its structure. In order to present a thorough investigation on this 
object, it is reasonable to write the current with isospin J = in a general form: 



which includes two adjustable parameters ci and C2- In case of ci = 0, Eq. (20) reduces to pure 
ss state which is free of instanton [l^ so we can deal it with the conventional QCD sum rule 
as in section II. The calculated mass is shown in figure 7 from which we can read the mass with 
pure ss structure: 




(20) 



1413 ~ 1430MeV, 



We can also write the current in a more complex form as: 




uu + dd) — 2ss 
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FIG. 8: Mass of /o with quark structure [{uu + dd) — 2ss]/\/6 as function of Borel parameter 
includes instanton contributions. 



such that there will be instanton effects in the corresponding sum rule. The calculated mass 
with this quark content is shown in figure 8 from which the following mass is read: 



m 



fo 



1433 ~ 1457MeV, 



In fact, based on the current given in Eq. (20), we can derive a sum rule involving complete 
instanton contributions induced by the quarks which depends on the two adjustable parameters 
Ci and C2. Unfortunately the results indicate the sum rule is not able to produce reasonable 
masses for the /o(1500) and /o(1710) by adjusting these two parameters — the masses are always 
much lower than the experimental ones. 

These simple calculations signal that from the viewpoint of pure quark content, one can 
not produce the /o(1500) in QCD sum rules even including complete instanton contributions, 
so we must look for other solution to this problem. As mentioned in Ref. j^, if we assume 
a qq structure, one concludes that /o(1500) is dominantly ss, while this assignment can not 
produce reasonable mass theoretically as we can see from previous paragraphs, but also leads to 
contradictions experimentally j^. There are some works 18 1 on this subject that take another 
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extreme: they try to produce /o(1500) under the assumption of a pure glueball content. But 
what is the reahstic structure of /o(1500) is still unknown. 

There is another viewpoint that the light nonet above IGeV can be identified as conventional 
qq states with some possible gluonic content, that is, there is mixing of the pure glueball with 
the nearby two = nfi and S = ss scalar mesons as first introduced in Ref. 



where 



nn = ^{uu+dd). Based on this model, Ref. 29| obtained the results that /o(1710) is dominated 



by ss content while there is roughly equal amounts of glue content in /o(1500). We have seen the 
key role of instanton in solving the puzzle on and oq, /o(1370), and explore this possibility 
in the assumed mixing of scalar meson and pure glueball in /o(1500). With this motivation, 
we modified the current of /o(1500) as mixing of quark and gluonic current:^ 

jrm, = Ass + BasG%G''^'', (21) 

and in this case the decay constant is defined as: 

{S\jmix\0) = mlfs- 

where A, B are both real, and one should notice that the parameter A has dimension one of 
mass which insures the right dimension in the current. The parameters A and B accompany 
the Wilson coefficients of operators ss and UsG'^^^.G'^^'^ respectively, and are therefore renormal- 
ization scale dependent. Here we fix the renormalization scale of A and B so that they just are 
numbers in the following consideration. After this modification, there will be new contributions 
stemming from the glueball asG'^^^^G"'^^ OPE, the glueball instanton and the mixing instanton 
contribution which will be presented below. 

With the perturbative corrections and including nonperturbative terms up to dimension 



^ The renormalization-group invariant gluonic current has been used because the subleading perturbative 
effects will be included in the correlation function. 
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eight, the OPE of gluonic current is 
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and 



(a) = (^7JabcG;,G::^G^^,) = (0.27GeV2)(a,G2), 
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are the dimension-6 and dimension-8 gluonic condensates, respectively. Because there is both 
quark and gluon current, we have to use the unsubtracted dispersive relation for the gluonic 
correlation function in order to be consistent with of the whole correlation function. Applying 
the dispersion relation, and after subtracting the continuum contribution and taking the Borel 



transform, the glueball contribution is obtained 
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The contribution of glueball instanton after subtracting continuum is given by 30l. l3l|: 

n^^'"^^*(so,M2) = -2Vr2p, / " dse-^s^J2{pvG)Y2{py^), 
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(24) 



17 



1.52 
1.50 



1.48 



1.46 



1.44 



1.42 



1.40 



2.00 2.05 



2.10 2.15 2.20 

M'(GeV') 



2.25 2.30 



FIG. 9: Mass of /o with structure O.OGeVss + a^G^^G"'^'^ as function of Borel parameter include 
glueball instanton and mixing instanton contributions. 



where J2 and Y2 are Bessel and Neumann functions, respectively. 

We have in dep endently verified the following instanton contribution to the mixed correlator 



321: 



923 rso 



(25) 



Now we have determined all the terms induced by the current given by Eq. (28). It is convenient 
to write the whole results in a compact form as follows: 



,4 ^2 „„„r "^Si 
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(26) 



X 



where X denotes 



X = \ {ss, OPE}, {GB, OPE}, {GB, inst}, {mix, inst} 



and we have absorbed the two parameters A and B in the 11^ 's for convenience. Taking the 
same algorithm as the previous section one can obtain immediately the mass corresponding to 
the current given in Eq. (21). Assigning A = O.QGeV and 5 = 1 in Eq. (21), corresponding 
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FIG. 10: Decay constant of /o with structure O.OGeVss + asC^^yG""^^ as function of Borel parameter 
include glueball instanton and mixing instanton contributions. 

to a large glueball content (since the energy scale is ~ IGeV), the calculated mass with this 
underlying structure is mj^ = 1492 ~ 1504MeV which is very close to the experimental result 
/o(1500). 

After obtaining the mass, we can deduce the decay constant from Eq. (26) 

//o(i50o)(lGeV) = 1.69GeV. 

There is a strong enhancement of the decay constant after involving glueball and related in- 
stanton contributions compared with other multiplets. A similarly enhancement of fc in pure 
glueball state including instanton effects was found in instanton vacuum model calculation 



. This strong enhancement also observed in the work of H. Forkel [30|] which gave a value 
fa = 1.14GeV when the glueball instanton contribution was included in the pure glueball cor- 
relation function using a unsubtracted dispersive relation. The results for our analysis of the 
mass and decay constant of this mixed state are shown in figure 9 and figure 10. 

Finally we turn to the last state /o(1710). Generally it is assumed this state dominated by 
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the ss content, so we can write the current as: 

J = A'ss + B'a^G^G''^'', (27) 
subjected the following orthogonality condition: 

(0|j|/o(1500)) = 0. 

This orthogonal condition is insignificant here since we are not able to get a value agreeable 
with /o(1710) whatever the values of A' and B' are chosen . This can be understood intuitively 
that the threshold sq = (1.9GeV)^, which is adopted here only for the states with the masses 
around 1450GeV, is too low to reproduce such a large mass. 

IV. CONCLUSIONS 

In this work we have studied the mass and decay constant of the light nonet ao, K^, and /o 
within the framework of QCD sum rule with and without instanton contributions. Our main 
results are as follows: 

1. In the conventional QCD sum rule, the masses of this nonet are degenerate, the calculated 
mass of Kq is larger than the oq for the same threshold and same Borel window. 

2. When we include instanton contributions in the sum rule, the masses of the nonet can be 
well separated, and the mass of Kq and ao agrees well with the observed results. The results 
suggest the underlying structure: ^0(1430) is sd, ao(1450) is -^{uu — dd), and /o(1370) is 
-^{uu + dd). For the /o(1500), our results suggest there is considerable glueball content in its 
underlying structure. The decay constant of /o(1500) enhanced considerably after this gluonic- 
content improvement compared with other multiplets. 

3. With a mixing current and the threshold and Borel window common to the multiplet, 
we cannot obtain the mass of /o(1710). One reason might be that the threshold suitable for 
JC(1430) is too low for /o(1710). 
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